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We theoretically study a doped graphene ribbon suspended over a trench and subject to an ac- 
electrical field polarized perpendicnlarly to the graphene plane. In such a system, the external ac- 
field is coupled to the relatively slow mechanical vibrations via plasmonic oscillations in the isolated 
graphene sheet. We show that the electrical field generates an effective pumping of the mechanical 
modes. It is demonstrated that in the case of underdamped plasma oscillation, a peculiar kind of 
geometrical resonance of the mechanical and plasma oscillations appear. Namely the efficiency of 
pumping significantly increases when the wave number of the mechanical mode is in close agreement 
with the wave number of the plasma waves. The intensity of the pumping increases with the wave 
number of the mode. This phenomenon allows selective actuation of different mechanical modes 
although the driving field is homogeneous. 

PACS numbers: 73.22.Pr, 81.07.0j,46.70.De 


Since the first graphene sample was isolated and stud¬ 
ied experimentally, the experimental and theoretical 
work on the 2D material has grown tremendously due 
to its many extraordinary propertiesJ^ The high mobil¬ 
ity, low mass and mechanical strength of graphene makes 
it well suited as the basis of nanoelectromechanical res¬ 
onators. The frequency tunability and high quality fac¬ 
tor of graphene based resonators make them promising 
for e.g., mass sensing^ and filtering applications.- To ac¬ 
tuate the nanomechanical resonators, different principal 
schemas are utilized. First of all, mechanical oscilla¬ 
tions can be initiated by applying an electrical field at 
resonance frequency with the mechanical vibration— ^2 
Another method which is utilized to control mechani¬ 
cal motion exploits the radiation pressure induced by an 
electromagnetic field in an optomechanical cavity— In 
this case the external frequency, at which the system is 
driven, is nonresonant with the relatively low mechani¬ 
cal frequency. Nonresonant excitation of the mechanical 
vibrations can also be achieved by integration of the me¬ 
chanical resonator in an electrical LC-circuit— S In both 
cases, the force acting on the mechanical subsystem is de¬ 
termined by the detuning of the external frequency and 
the resonance frequency of the cavity or external LC- 
circuit. The resonance frequency depends on the mechan¬ 
ical displacement which induces an electro-mechanical 
time-delayed backaction. The backaction generates an 
effective pumping (or damping) of the mechanical vibra¬ 
tions. Therefore, it is possible both to excite and cool the 
resonator. These phenomena have been demonstrated for 
many system s and for graphene based resonators in 
particular— Recently it was shown that similar effects 
can be achieved by integrating the resonator into an RC- 
circuit— 2 In this description, the actuation mechanism 
was due to the time-delayed overdamped charge response 
rather than coupling via a resonant high-frequency mode. 

In this article, we show that nonresonant excitation of 
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FIG. 1: An isolated graphene membrane is suspended over a 
cavity inside a wave-guide. The suspended part of the mem¬ 
brane is free to perform vertical vibrations. A electrostatic 
ac-field is applied polarized perpendicularly to the graphene 
sheet. 


mechanical vibrational modes can be achieved also for an 
isolated graphene membrane via its internal charge dy¬ 
namics. We will demonstrate that the nonresonant actu¬ 
ation mechanism presented here enables selective actuate 
of different mechanical modes, even antisymmetric ones. 
The intensity of actuation increases with mode number in 
contrast to the optomemechanical and electrical pump¬ 
ing mechanisms mentioned above where predominantly 
actuation of the fundamental mode takes place. 


I. MODEL 

The sketch of the system under consideration is pre¬ 
sented in Fig. [U and comprises the isolated graphene 
ribbon suspended over a trench with width 1. The sys¬ 
tem is placed in a wave guide. Electromagnetic waves, 
with wave length much larger than I travel inside the 
wave guide along the x axis. The field is assumed to be 
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homogeneous in the trench and screening effects from the 
wave guide is neglected. The wave is polarized perpen¬ 
dicularly to the flat membrane and induces an electrical 
field along the membrane only when it is deflected from 
its flat position. The induced field generate electronic 
charge waves in the graphene sheet. Simultaneously, the 
electrical field exerts a force on the suspended part when 
it is charged and provides a feedback coupling between 
the electronic and the mechanical subsystems. 

To analyse this feedback we model the free vibrating 
part of the ribbon as an elastic membrane. For sim¬ 
plicity we consider the membrane to be infinite in the 
y-direction. In this limit, we assume that the membrane 
deflection U{x,y,t) = U{x,t), charge density Q(x,y,t) = 
g{x,t) and current density jx{x,y,t) = j{x,t) are uni¬ 
form along the trench. We disregard the geometric non¬ 
linearity of the graphene membrane since it does not af¬ 
fect the nonresonant phenomenon discussed in this paper 
and can be neglected at small amplitude of oscillation 

C/<h 

Under these assumptions the dynamical equations for 
the flexural out-of-plane modes become 


f d To K d^\ , 

V dt dx^ ^ dx* ) 


g{x,t) 


m 

( 1 ) 


with intrinsic mechanical damping 7 , built-in tensile 
stress To, bending rigidity of graphene k, electrical field 
in the wave guide E{t) = Tq cos(nt) and 2D-mass den¬ 
sity of graphene pm- The corresponding boundary condi¬ 
tions of the clamping are U{x,t) = 0 and U'{x,t) = 0 at 
X = ±Z/2. The membrane deflection U{x,t) can be pre¬ 
sented as a superposition of the vibrational eigenmodes 


U(x,t) =^Un(t)fn{i) (2) 

n—1 

with dimensionless spatial coordinate ^ = x/l, stretching¬ 
bending ratio = k/tt'^PTo and fn{0 is the normalized 
spatial profile of the flexural eiqenmodes. 

To describe the charge dynamics of the electronic sub¬ 
system we will use a simple hydrodynamic approach d^d^ 
We will consider monopolar electronic plasma where the 
Fermi energy Ep is much greater than temperature and 
ftU. Within this approach the charge evolution is de¬ 
scribed by 

^g{x,t) =-^j{x,t), (3) 

d 1 

—j{x,t)+pj{x,t) = -E{t,x), (4) 

where 1/C = e^Ep/b?Tr and v is the scattering frequency. 
The electrical field along the ribbon consists of one ex¬ 
ternal and one internal contribution 

E{t,x) = E{t)^U{x,t) + (5) 

ox ZTTeo J X — xi 


TABLE I: Dispersion relation for plasma and mechanical vi¬ 
brations. The continuons wave number q corresponds to a 
plasma wave length A = 2Z/|g|. 


Electronic 

Mechanical 

uj{k) = lXp^J\q\ 

, . _ e Ep 
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The first term in Eq. ([S]) describes the external electrical 
field induced along the membrane when it deflects from 
its flat position. The second term describes the internal 
non-local electrostatic field due to charge redistribution. 

The time-scales of the system is obtained by consider a 
typical experimental situation where we take I ~ 10 /im, 
Ep Ri 1 meV, end Tq ~ O.lN/m which gives b ~ 10“^. 
Under such conditions the characteristic mechanical fre¬ 
quency ujm ~ 100 MHz and the characteristic plasma fre¬ 
quency ujp u PS I THz are well separated. Further, we 
consider high-frequency external driving H ^ Wp ^ wm- 
The strength of the electromechanical coupling generated 
by the external field is characterized by the coupling fre¬ 
quency top = Eo\/eo/2TTlpm- We will consider low am¬ 
plitude external field so that uip is the smallest frequency 
iOp « OJm << ^ 2 . 


II. EFFECTIVE MECHANICAL DYNAMICS 


To get the coupled dynamics for the amplitudes u„(t) 
and charge density g{^,t) we combine Eqs.dS]) and (|4]) 
and obtain 


Unit) + JUnit) + ulfKluit) = ^^(/n(C), £'(C, 0) 


d r 

gii,t)+vgii,t) - J 


gi^i,t)d^i 


Un{t)-^fniO, 
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where (fniO^ t)) denotes projection of the charge dis¬ 
tribution on the spatial mode function /„. The charac¬ 
teristic mechanical and plasma frequencies ujm and tOp 
are defined in table HI 

The electrostatic forces acting on the vibrational 
modes can be expressed by substituting an integral ex¬ 
pression for the charge density described by Eq. © 
in the right hand side of Eq. (|6l). The forces En{t) = 
ddit){fniO> Qiift))/pm can then be formulated as 

Tn it)= ( 8 ) 

Lm / dtiGnmit - h) COsinit - ti))Umitl) + 
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with 


Gnm (^) — 


2 e sin (^^uip\q\ — 

Y^w2|g| - 1/2 


w„ 


i(9)dg. 


(9) 

where w„m(g) = -(7rg)^(/„(0> e“^‘^)(e “^‘^,/m(C))- 

The electrostatic forces on the form Eq. m introduce 
linear feedback on the mechanical motion. The feedback 
on mode n is direct back to itself via G„„(t) but the 
feedback also couples different modes via Gnm{t) n ^ m. 
We want to note that the subsets of odd and even modes 
do not couple. However, since the coupling strength e = 
We/wm is assumed to be small we disregard the coupling 
between modes since it will affect the mechanics only to 
fourth order in e. The system of equations Eq. (jS]) then 
decouples to independent single mode oscillators. 

The dynamics is further simplified since we consider 
the high-frequency regime of the driving frequency H ^ 
uip >> wmjWe- As we will see later, under such condi¬ 
tions only modes with Kn ~ (Q./ujp)'^ play an important 
role in the membrane dynamics. We seek the time evolu¬ 
tion of the amplitudes Unit) in the form of perturbation 
series 


= 


E 2 m amnt 

c D ^n,m 

m —0 
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here e = ojm/^ << 1 and Un,mit) are slow on the time 
scale The second term in the right hand side of Eq. 
(|5|) gives corrections of the order e^e^. We will neglect 
corrections of this order of smallness and take Unit) ~ 
Un,oit)- In these approximations, the dynamics of mode 
n is governed by 

Unit) + 'yilnit) + uljKlunit) = (11) 

/ dfi G„n(t - fi) cos(H(t - fi))u„(ti). (12) 

J —OO 
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FIG. 2: Spatial geometrical resonance between the charge 
oscillations and the mechanical mode functions occur in the 
vicinity of g = Kn = n. Wave number Kn = n corresponds 
to a mechanical mode with wave length X — 21/n. 


here Q = ujm/ l is the characteristic Q-factor and we have 
introduced the dimensionless frequencies Cjm = 
and XI = Xl/ujp. The functions AiXl,q) and r]iXl,q) are 
given by 


Ain,q) 


- . 




a) - ) 


(16) 

(17) 


and Wnniq) = ( 7 rq)^|(e™^'?,/„(^))p. To calculate Wniq) 
we neglect the bending rigidity of the graphene in com¬ 
parison with the clamping tension since 6 << 1. Un¬ 
der this condition we can take fniO — sin( 7 rn^ -|- 
TTn/2)9il/2 — \^\)/^/2 with Kn = n. These forms for 
the mode shapes give 


Wnniq) = 


2 ( 7 ^n^ sin'^( 7 r(|q| -|-n)/ 2 ) 
(<72 - n2)4 ■ 


(18) 


III. SELECTIVE MODE ACTUATION 


The dispersion relation which characterises the time 
evolution of the n:th mode can be obtained by the Ansatz 
Unit) = exp(ia;„f). Substituting this form in Eq. (fTTl) we 
obtain 


~ Wn T -|- Uj\[Kn — Uj'^GnniWn', Xt) (13) 

G„Ti(a;;0) = — / dt exp(—fwt)G„„(t) cos(ilt) (14) 

J — OO 

Solving Eq. m we arrive at the following approximation 
(with an accuracy of e^) for the complex frequencies a;„ 


w„ — Wm 


Kn + |-A„(11) + ^ (^ + 


A„(0) = J Ain,q)wnniq)dq, 
r]ni^) = / r]iXl,q)wnniq)dq. 


(15) 


The shift of the complex frequencies occure at strong 
coupling between the mechanical modes and the charge 
waves. The charge waves are mainly generated from the 
regions close to the clamping since the gradient of the 
electrical field along the membrane is biggest in these ar¬ 
eas. The coupling is strong when the wave length of the 
generated charge waves is in close agreement with the me¬ 
chanical wave length. This can be seen in Fig. [2]since the 
functions Wniq) has a sharp maximum in the vicinity of 
Igl ^ n. Therefore we have a peculiar kind of spatial geo¬ 
metrical resonance but a nonresonant phenomenon in the 
time domain. Simultaneously, at small plasma damping 
E << 1, the functions r]iCt,q) and Airt,q) dramatically 
increase when jgl = Xl^. 

The shift of the mechanical damping is qualita¬ 
tively described by the normalized damping coefficient 
77 „(H)/n^, Fig. 131 The damping coefficient of the nth 
mechanical mode rjniXl/ujp) becomes negative at H = 
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FIG. 3: The damping coefficient for damping ratio v equal 
3, 1 and 1/3 in (a), (b) and (c), respectively. A region of 
negative damping coefficient is present at driving frequencies 
Q. > uip. The mechanical modes can selectively be driven if 
f' < 1 due to the geometric resonance of plasma and mechan¬ 
ical oscillations. 


this, the phenomenon is nonresonant in the time domain. 

A vibrational mode will become mechanically actuated 
if the effective pumping generated by the high-frequency 
external field overcomes the intrinsic mechanical damp¬ 
ing of the mode. To actuate a mode the driving fre¬ 
quency has to be in the region where the electromechan¬ 
ical cpupling gives negative damping and the amplitude 
of the external field has to exceed a critical value E^{n). 
Above the critical value the mechanical vibration is un¬ 
stable and will be saturated by nonlinear effects. The 
field strength needed to achieve this can be estimated by 
using Pm Ri 0.7 mg/m^, assuming the quality factor to be 
Q = 10® and damping ratio ly/ujp = 1/3. This gives an 
estimate of the critical field strength for the fundamental 
mode El r; 30 V//rm, at the optimal driving frequency 

Selective actuation of vibrational modes is possible 
when the overlap of the minima peaks is small, Fig. [3]:. 
This possibility is remarkable since the applied electri¬ 
cal field is homogeneous. It is interesting to note that in 
contrast to optomechanical excitation where only sym¬ 
metric modes are actuated and the strength of pumping 
decreases with mode number, whereas in our system also 
antisymmetric modes can be actuated and the strength 
of pumping increases with mode number. 

From the above analysis it follows that pronounced se¬ 
lective nonresonant excitation of the mechanical modes is 
achievable for z/ <C 1. However, there are natural restric¬ 
tions for V and n which come from the applicability of the 
hydrodynamic description of the charge dynamics used in 
this paper. The hydrodynamic equations are not valid in 
the ballistic regime of electronic propagation. To analyze 
the range of parameters where our approach is valid, it is 
convenient to introduce the effective electron mean free 
path Isc = vfIv. We then have the following expression 
for the damping ratio v = vjuip = 0.271 j(Iscy/kpl) — 
0.071/Isc- The hydrodynamic approach fails when the 
electron mean free path exceeds characteristic space vari¬ 
ations in the system Igc > l/n. Therefore, the model 
formulated here is not valid for n > 14z>. 


(f') ^ ujp and reaches its minima rj™'^{y) at the min¬ 
ima frequency If the plasma oscillation is over¬ 

damped Fig. [3^, the characteristic width of the minima 
is much greater than the distance between the minima 
frequencies while ??J/'*"'(3) r; —ri^. In the under¬ 

damped situation Fig. [St, the distances between the min¬ 
ima frequencies become greater than the width of the 
minima and pumping strength 77j/’™(l/3 ) r; —lOn^. It 
should be particulary emphasised that the distance be¬ 
tween minima frequency as well as minima widths are 
three order of magnitude greater than, and independent 
of, the characteristic mechanical frequency. Because of 


IV. CONCLUSIONS 

To conclude, we have shown that the internal charge 
dynamics in a suspended isolated graphene sheet can be 
utilized to selectively actuate vibrational modes by a non¬ 
resonant homogeneous external field. The phenomenon 
occurs when the external field induces plasma oscillations 
with a wave length comparable to the wave length of the 
spatial profile of the vibrational mode. Different modes 
can then be selectively driven via this geometrical reso¬ 
nance, if the plasma oscillations are underdamped. 
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